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Abstract 



We study nuclear reactions induced by virtual photons associated with 
Lorentz-boosted Coulomb fields of ultrarelativistic heavy ions. Evapora- 
tion, fission and multifragmentation mechanisms are included in a new 
RELDIS code, which describes the deexcitation of residual nuclei formed 
after single and double photon absorption in peripheral heavy-ion colli- 
sions. Partial cross sections for different dissociation channels, including 
the multiple neutron emission ones, are calculated and compared with 
data when available. Rapidity and transverse momentum distributions of 
nucleons, nuclear fragments and pions, produced electromagnetically, are 
also calculated. These results provide important information for designing 
large-rapidity detectors and zero-degree calorimeters at RHIC and LHC. 
The electromagnetic dissociation of nuclei imposes some constrains on the 
investigation of exotic particle production in 77 fusion reactions. 
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1 Introduction 



There exist many mechanisms of particle production, nucleon and fragment emis- 
sion in nuclear collisions at ultrarelativistic energies. They can be divided into 
the following ample categories. The first category contains collisions with di- 
rect nuclear overlap. In addition to nucleon-nucleon or parton-parton collisions, 
which are the most obvious sources of secondary hadrons, collective mechanisms 
were considered for baryon-antibaryon pair |l] or dilepton || production. The 
second category deals with "distant" electromagnetic collisions without nuclear 
overlap. In such a situation the interaction can be treated in terms of virtual 
photons emitted by one nucleus and absorbed by another collision partner. A 
subsequent photonuclear (7 A) reaction leads to a so-called Electromagnetic Dis- 
sociation (ED) of nuclei || |J. The virtual photon flux and consequently the rate 
of such 7 A events is proportional to Z 2 , where Z is the charge of the first nucleus. 
The ED process is widely used for studying the properties of Giant Resonances 
(GR) in nuclei, including multiphonon excitations ||. 

In Ref. || the description of nuclear photoabsorption was extended to the 
photon energies much above the GR region, where the exitation of individual 
nucleons and multiple pion production become important. The model of elec- 
tromagnetic dissociation processes taking into account these high-energy photon 
absorption channels was constructed in Ref. 0. According to this model, the 
fast hadrons produced after the photon absorption initiate a cascade of subse- 
quent collisions with the intranuclear nucleons leading to the heating of a residual 
nucleus. At a later stage the nucleus undergoes de-excitation by means of the 
evaporation of nucleons and lightest fragments, binary fission or multifragmen- 
tation. The latter process takes place when the excitation energy is comparable 
with the total binding energy of a residual nucleus f?j. 

Another interesting mechanism in the second category is the photon-photon 
collision or 77 fusion leading to the particle production at mid-rapidity [§, |9], [TtJ. 
The rate of these processes is proportional to Z A . The effective 77 luminosity 
at the LHC collider is expected to be comparable with that at a possible future 
Photon Linear Collider [|TT|j . 

It is worthwhile to mention that the processes of the second category have 
attracted a special attention in the last few years. This interest is generated by 
the construction of the new heavy-ion colliders: RHIC at Brookhaven (Au + Au 
with v/i = 200 GeV) and LHC at CERN (Pb + Pb with ^ = 5500 GeV). It is 
expected that the ED processes will play an important role for these energies and 
especially for the ions with high Z, due to the above-mentioned Z 2 -dependence. 
The predicted cross section for the electromagnetic dissociation at LHC, oed, is 
so large (~ 200 b) that together with the electron capture process (~ 100 b), 
it will reduce significantly the life time of relativistic heavy-ion beams |12[ as 
compared with the proton ones. The ED reaction rate at LHC is expected to be 
30 times higher than that for the nuclear interactions (a nuc « 7 b). But of course 
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the ED reactions are much less violent then the nuclear ones. 

Hadron production in photon-photon and photon-pomeron collisions, as well 



as in double-pomeron exchange were studied in Ref . |T3| . All these reactions are 
characterised by a large rapidity gap separating the remnants of the colliding 
hadrons or heavy ions from the particles produced in the central rapidity region. 
High-mass diffraction reactions were simulated by means of the Monte-Carlo 
event generator Phojet On the other hand, there are no "ready-to-use" event 
generators for simulating relatively "soft" 'jA processes, where nuclear degrees of 
freedom play an important role. In particular, such simulations are needed for 
estimating the probabilities of different disintegration channels of the colliding 
nuclei. The properties of the ED processes have to be well understood in order 
to distinguish them from very peripheral nuclear collisions. 

In the present paper we consider the whole set of mechanisms of Coulomb 
excitation and subsequent dissociation of relativistic heavy-ions (RHI). In Sec. 2 
the model is presented and compared with existing experimental data on few 
neutron removal from Au and Pb nuclei by virtual and real photons. In Sec. 3 we 
discuss the ED processes which are expected to be important at the future heavy- 
ion colliders, RHIC and LHC. We present the model predictions for multiplicity 
and transverse momentum distributions of neutrons. These results are important 
for designing zero-degree calorimeters for beam luminosity monitoring. Rapidity 
and transverse momentum distributions of ir + and tt~ mesons, fission products 
and intermediate mass fragments are also discussed in Sec. 3. Sec. 4 is reserved for 
discussion and conclusions. 



2 Model of electromagnetic dissociation of heavy 
nuclei 

2.1 Single- and double-photon absorption 

Let us consider the Coulomb excitation of an ultrarelativistic projectile. Accord- 
ing to the Weizsacker- Williams (WW) method the spectrum of virtual photons 
from a stationary target of charge Z t as seen by a projectile moving with velocity 
(3 = v/c at impact parameter b is expressed as |15| : 

N ^ h ) = ^^(^« + ^o(x)), (1) 

where a is the fine structure constant, K and Ki are the modified Bessel func- 
tions of zero and first order, x = E^ 1 b/( , yf3hc). It is worth noting that the distri- 
bution (0) is written in the projectile rest frame. This means that for colliding 
beams one should use the Lorentz factor of one collision partner in the rest frame 
of another partner, 7 = 2^ eam — 1, where 7f> eam is the Lorentz factor of each 
beam. 
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The projectile can be exited by absorbing one or more virtual photons. Mul- 
tiple excitations of the giant dipole resonance (GDR) up to the fourth order were 



considered in Refs. EL 16|. Since the nth-order cross section behaves approxi- 



mately as («Zj 2 )", one should expect the maximum effect for the excitation of 



J t 

238yj projectile in the field of 238 U target at ultrarelativistic energies. As was 
found in Ref . , even in this case the cross sections of the third and fourth order 
processes are quite small in comparison with that for the second order excitation: 
cr (3)y (7 (2) ^ o.ll and cr^/en 2 ) m 0.01. Thus we restrict our consideration by the 
first and second order processes. 

The mean number of photons absorbed by the projectile of mass A p in a 
collision at impact parameter b is defined by 

oo 

m(b)= J iV(£ 7 ,6)a Ap (£ 7 )d£ 7 , (2) 

where cta p (E^) is the appropriate photoabsorption cross section, either measured 
for this projectile nucleus with real photons or calculated within a model. For 
the cases of interest, i.e. for Au and Pb nuclei, we use the data set published 



in Refs. fL7\, ffq, 19, 20], including a most recent one [21| where corresponding 



photoabsorption cross sections were measured. In the GDR region one can use 



a set of Lorentz-line approximations given in the review [22|. According to pre 



scriptions of Ref. p0 |, the photoabsorption cross sections obtained in Refs. |~7| 



and |n| for Au and Pb nuclei are rescaled by factors 0.93 and 1.22, respectively. 

Following Llope and Braun-Munzinger |16[ we assume that at an impact pa- 
rameter b the probability of the multi-photon absorption is given by the Poisson 
distribution with the mean multiplicity m(b) defined by Eq. ([|). Then the proba- 
bility for absorbing exactly one photon of any energy during a collision at impact 
parameter b is equal to: 

pW(fo) = m {b)e- m{b \ (3) 
and the probability for absorbing exactly two photons is: 

pi 2 ){b) = r^) e -m ( 6)_ (4) 

The normalized probability density that in a collision at impact parameter b 
the absorbed photon has an energy E\ is given by: 

m(b) 

Analogously, the probability density that in a second order process two photons 
have energies E\ and E 2 is: 

(aj/p P w NjE^A^NjE^aA^E,) 

<T '{Ei, E 2, b) = 277a • 6 

m z {b) 
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To consider particular ED channels i.e. neutron emission, fission, multifrag- 
mentation or pion production, which are studied in the present paper, we define 
the differential cross sections for the first and second order processes as: 



da 



(i) 



DC 



and 



dE l 

b 



dV (2) 



2tt J bdbP^(b)q^ 1 \E 1 ,b)ft\E 1 ) (7) 



dE 1 dE 2 



= 2n J bdbP^\b)q^{E u E 2 ,b)ff\E l ,E 2 ) (8) 

Here // (.Ex) and (Ei, E 2 ) are the branching ratios for the considered chan- 
nel i. 

Finally, the integral cross sections are calculated as: 

oo 

a^ = J dE 1 N( 1 \E 1 )a Ap (E 1 )f^(E 1 ) 1 (9) 



(2) 

07 



J J dE 1 dE 2 N^(E 1 ,E 2 )a Ap (E 1 )a Ap (E 2 )fP(E l ,E 2 ). (10) 

Here we introduce the single and double photon spectral functions which appear 
after the integration over impact parameters; 

oo 

N w (E 1 ) = 2n J bdbe- m{b) N(E u b), (11) 

OO 

N {2 \E U E 2 ) = 7T J bdbe- m ^N(E u b)N(E 2 ,b). (12) 

In the above expressions b m i n is the minimal value of the impact parameter 
which corresponds to the onset of nuclear interaction. The integration over the 
photon energy starts from E m i n ~ 7 MeV which corresponds to the threshold of 
nuclear dissociation in photonuclear reactions. Due to the long-range nature of 
the electromagnetic forces, numerical calculations of the integrals ([TT|) and (|i~2|) 
should be undertaken with care. A reasonable accuracy can be obtained if one 
splits the whole region of integration into two intervals: [b min , b cut ] and [b cut , oo). 
The value of b cut is determined by the condition e~ m ^ bcut ^ ~ 1 which allows one 
to omit the exponential factor in the second interval. After that the first integral 
in Eq. (|IT]) is calculated numerically, while the second one can be evaluated 
analytically, leading to a well-known expression |Tj|. The same splitting can be 
made in the integral (O). Using the asymptotic behaviour at x 3> 1: 

tfi(x) « v^72xe- x , (13) 
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one obtains the exponential integral which can be found in mathematical tables. 
With the described method one can avoid problems with too many nodes in the 
numerical integration over the wide range of impact parameters. 

The values of and for branching ratios of partial channels of photonu- 
clear dissociation induced by real and virtual photons were calculated within the 
model of Ref. |J which takes into account the fast cascade stage of the reaction 
as well as evaporation, fission and multifragmentation processes for a residual 
nucleus. 

2.2 Excitation and decay of residual nuclei 

Relative probabilities of the de-excitation processes are determined by the exci- 
tation energy E* of a residual nucleus formed after the completion of the fast 
cascade stage of the photonuclear reaction. In the following we consider step-by- 
step all these mechanisms with a special attention to the amount of energy which 
is transformed on average into the internal excitation of the nucleus. 

Depending on the virtual photon energy different processes may contribute 
to the energy deposition. When a nucleus absorbs one or two virtual photons 
in the GR region, 6 < E 1 < 30 MeV, their energies are completely transformed 
into the excitation energy E*. For preactinide nuclei like Au and Pb, whose 
fission thresholds are about 30 MeV, the de-excitation proceeds mainly through 
the evaporation of neutrons, since their separation energies are only around 7 
MeV. Due to a high Coulomb barrier in heavy nuclei, the proton emission is 
suppressed in the GR region. A detailed experimental investigation of (7, n), 
(7, 2n), (7, 3n), and (7, 4n) reactions on 197 Au and 208 Pb nuclei in the GR region 
was performed in Ref. [JI7| . The presence of the reaction (7,4n) means that a 
nucleus has received E* up to 30 MeV. For simulating the particle emission in 
the evaporation process we used the standard Weisskopf evaporation scheme J7| 
with the assumption E* = E 1 . 

Starting from = 30 MeV, where the quasideuteron mechanism becomes 
important and up to the single pion production threshold at E n = 140 MeV, only 
a part of the photon energy is converted on average into the excitation energy 
E* of the compound nucleus. The rest of the photon energy is taken away by 
the fast nucleons originating from the absorbing pair. The energy deposition is 
simulated by the IntraNuclear Cascade (INC) model (see Refs. 0, |6[ for details). 
The two-nucleon absorption cross section of a photon on a heavy nucleus, a a, is 
taken from the quasideuteron model of Ref. |18[ : 

a Q A D = kZ(l - Z/A)a e d xch - (14) 

Here a^ xch is the meson exchange part of the cross section for the deuteron pho- 
todisintegration, / yd — > np. This cross section was calculated in Ref. ||23|| . In 
Eq. (|Hp A and Z are the mass and charge numbers of the relevant nucleus 
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and k ~11 [18|] is an empirical constant. The angular distribution of nucleons 
in the reaction 7c? — > np was approximated in accordance with empirical fits 
from Ref. fl24|. Although the cross section ad decreases strongly with the pho- 



ton energy, the two-nucleon absorption mechanism competes noticeably with the 
single- nucleon photoabsorption channel up to E 1 ~0.5 GeV. 

As found in Ref. |18| , up to 12 neutrons can be emitted by photo-excited 



Pb nucleus. Assuming that the emission comes from the compound nucleus one 
can conclude that (E*) < 70 — 80 MeV. The calculations of Ref. |2(J for pho- 
toabsorption on 197 Au and 208 Pb nuclei, performed in the framework of another 



version of INC model [p4], demonstrate that (E*) does not exceed 80 MeV for 
30 < £ 7 < 140 MeV. In this £ 7 region the fissility P f of 197 Au and 208 Pb nuclei 
becomes noticeable (Pf ~ 0.01 and 0.1, respectively), so one should expect the 
evaporation-fission competition. 

Above the single pion production threshold at E 1 = 140 MeV the A-isobar 
excitation of individual nucleons dominates. An interesting effect was noticed in 
Ref. j2~B|. Near the threshold a slow pion has a small interaction cross section 



with nucleons and therefore, has a high probability to carry away a large part 
(m m n ) of the photon energy. Only at E 1 ~ 200 MeV the average value (E*) 
starts to increase again. During the fast cascade stage the initial nucleus looses on 
average 1-2 nucleons and acquires the excitation energy (E*) ~ 100 MeV, which 
is sufficient for evaporation of many neutrons or fission. The evaporation may 
take place before of after the fission. In some cases fission fragments evaporate 
neutrons as well. 

Above the two-pion production threshold at E^ ~ 400 MeV the photon- 
nucleon (7AQ interaction becomes more complicated because of many possible 
final states. Below we use a phenomenological model for exclusive description of 
the 7iV-interaction developed in Ref. f| . The model includes both the resonance 
contribution from the two-body channels, 7iV — > nB* and 7iV — > M*N (B* and 
M* being baryon and meson resonances), and the non- resonant statistical con- 
tribution from the multi-body channels 7iV — > inN (2 < i < 8). Subsequently, 
the fast hadrons produced in a primary 7iV-interaction initiate a cascade of suc- 
cessive hadron-nucleon collisions inside the nucleus. Finally, when the photon 
energy reaches the value of several GeV, the multiple pion production becomes 
the main process. A large number (~ 80) of many-body subchannels is included 
in the calculation. 

Fig. m shows the fraction of E 1 that on average is transformed into E* as 
well as average values of E* per nucleon of residual nucleus. As it was deduced 
from the experimental data [^5 |, E* = A3 A ± 5 MeV for photoabsorption on 
lead at E 1 = 70 MeV. This value agrees well with the prediction of our model 
(see Fig. |TJa) . In Fig. |l|a one can see the above-mentioned decrease of (E*)/Ery 
above E 1 = 140 MeV. As seen from Fig. when the virtual photon energy 
grows from the GR region to several GeV, the nature of the Coulomb excitation 
changes from the collective nuclear degrees of freedom to the excitation of a single 
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nucleoli inside the nucleus. In the latter case up to 95% of the photon energy 
is released in the form of fast particles leaving the nucleus. Nevertheless, the 
remaining energy deposited in the residual nucleus is sufficient for evaporating 
many neutrons. The moments of the neutron multiplicity distribution are shown 
in Fig. [I], parts c and d. Calculated moments are in a satisfactory agreement 
with the experimental values for the real photon absorption on Pb nucleus |p5 |. 



Below we use the same model for the case of the virtual photon absorption in 
electromagnetic interactions of ultrarelativistic heavy ions. 

It should be stressed however, that Fig. [I] shows only the average values, while 
a part of the interactions initiated by "hard" photons will lead to a stronger 
heating of the nucleus by a multi-pion system. This phenomenon was studied in 
detail in Ref. [H] were it was shown that the yield of events with high excitation 
energy, resulting in the multifragment break-up of residual nuclei, might be quite 
sizable, up to 10 -r- 15% of the total ED cross section. These events can be 
easily observed with properly designed forward-angle detectors at RHIC and LHC 
colliders. 

To include all the processes described above we have developed a specialized 
computer code RELDIS aimed at the Monte Carlo simulation of the Relativistic 
ELectromagnetic Dissociation of nuclei. The simulation begins with generating 
the single or double-photon absorption process. Then the intranuclear cascade 
model is used to calculate the fast particle emission and the characteristics of 
residual nuclei. Finally the de-excitation of thermalized residual nuclei is simu- 
lated according to the statistical model 0. 

2.3 Analysis of nucleon emission data 

Let us consider now how the calculations based on Weizsacker- Williams (WW) 
method describe the measured |27]] cross sections of In, 2n and 3n removal from 
197 Au nuclei excited electromagnetically. The branching ratios, and f\ 2 \ 
for partial channels of a photonuclear reaction are important ingredients of such 
calculations. The information concerning these ratios can be obtained from ex- 



periments with real photons [jD], |20|, 23[] or from models of photonuclear reactions 
analogous to the one used in the present paper. We review in brief the advantages 
and disadvantages of both these approaches. 



The calculations [27] within the framework of WW method, which include 
both isoscalar and isovector quadrupole excitations as well as multiphonon states, 
describe the one-neutron emission data fairly well. Using the harmonic-oscillator 
model one can derive a Poisson distribution for the excitation probabilities of 
multiphonon states. On the contrary, even the theory with these additional 
modifications is not able to describe the 2n removal cross section. The measured 
cross section is twice as large as the calculated value, which is determined mainly 
by the probability of the double-phonon excitation. Further modifications of the 
WW method, using a "soft-spheres" model and including nuclear-plus-Coulomb 



S 



processes 
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do not remove the discrepancy. The description of the reaction 
mechanism within the coupled-channel approach |29|] with proper account for the 
nuclear interaction is not able to improve the situation. Only the inclusion of 
anharmonicities in the internal Hamiltonian and nonlinear terms in the external 
field, proposed in Ref. JJUJ and later supported by other authors [[31], help to solve 
the problem. As shown in Ref. |30], [H]] , a small deviation from the harmonic model 
may lead to a doubling of the two-phonon excitation probability and to a good 
description of the data. 

The calculations of Refs. 



30, 31 



are based on the one- and two-neutron 
emission cross sections measured in experiments with real photons. Thus only a 
limited number of channels of electromagnetic dissociation was considered. The 
differences in the total GDR cross section measured by different groups are found 
to be on the level of 20% 



20 



There exists also an opinion [32] that for the 2n 
photoemission measurements the quality of the data is even much worse. Different 
interpretations of the photonuclear data lead to difficulties in determining the 
contribution of the double-photon absorption. Indeed, the double GDR state 
may decay via 2n and 3n channels and therefore its admixture depends on the 
corresponding partial cross sections. 

In the present paper we analyse the ED data [27] within the model of pho- 
tonuclear reactions introduced above. We do not use experimental data as an 
input for calculations, except of the information concerning the non-statistical 
contribution in the In emission channel. It should be reminded, that the statis- 
tical decay of the compound nucleus is not the only mechanism responsible for 
the decay of the collective GDR state. This state can also decay directly into 
a free nucleon and a hole in the residual nucleus, due to the two-body residual 
interaction in the final state (see Refs. [|T7|, Q for details). Such a direct neu- 
tron emission leaves the residual nucleus in a low-lying energy state where only a 
subsequent emission of photons is possible. This makes the emission of a second 
neutron impossible, leading to a suppression of the (7, 2n) channel in compar- 
ison with the pure statistical decay. We include this additional non-statistical 
contribution in the (7, n) channel on the basis of the experimental information 
concerning this decay mode. The total fraction of the non-statistical contribution 
was evaluated to about 35% and 20% for Au and Pb nuclei, respectively |17]]. We 
use these values in our simulations. The emission angles of non-statistical fast 
photoneutrons were generated according to the angular distributions measured 
in Ref. || . 

In our model of photonuclear reactions, implemented in the RELDIS code, 
we consider all possible decay modes, including binary fission and multifragmen- 
tation. It gives the yields and spectra of different particles in the ED reactions 
induced by real and virtual photons. Before discussing the predictions in the 
remaining part of the paper, let us consider first how the model describes the 
existing data on a few neutron removal from Au nuclei [27] . 

The results for the few-nucleon emission are summarized in Table |IJ One can 
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see that the model describes the 3n emission quite well but the agreement with 
the measured one- and two-neutron emission cross sections is not as good. At 
the same time the total neutron (n — 3n) emission cross section agrees with the 
experimental value. Due to the statistical nature of our model one can expect that 
it is not so accurate for the particle emission from collective states like GDR, but 
works better for the multiple emission from the compound nucleus. This tendency 
is clearly seen in Table [l]. The higher is the energy of colliding ions the greater is 
the number of emitted particles. Therefore, the quality of description of the low 
energy data p7| may be considered as a low limit for the accuracy of the model 
which is expected to work better with increasing energy. It should be noticed 
that a simple Poisson distributions for multiple excitations was adopted in our 
calculations. And even within this assumption the agreement with data pTJ is 



acceptable, in contrast with observations of Refs. |2j, p8[ . This discrepancy is 
due to slightly different values of the partial photonuclear cross sections (2n and 
3n) used in our calculations as compared with the papers pTj ^8 . 



The importance of the double-photon absorption mechanism is illustrated 
in Tab. [I] too. Although the total cross section of double-photon absorption, 
<j( 2 ) = 226 mb, is small compared to the single photon one, = 3546 mb, 
this contribution is very significant for the rare channels like (7,4n), (7,5n), 
(7,6n) and, especially, for (7,p3n), (7,p4n) and (7,p5n) channels. These listed 
channels are extremely sensitive to the presence of the double photon absorption 
mechanism. 



3 Predictions for ED at RHIC and LHC 

3.1 Total and partial ED cross sections 

Let us consider now the model predictions for the electromagnetic dissociation of 
ultrarelativistic Pb and Au nuclei on fixed targets at SPS and in colliding beams 
at RHIC and LHC. The total ED cross section, oed — + is obtained 
from Eqs. @ and ([U]) by taking = 1 and f\ 2 ' = 1. The Oed values for most 
interesting reactions are presented in Tab. |2|. 

The inclusive (multiplicity weighted) cross sections for emission of nucleons, 
nuclear fragments and pions in the electromagnetic dissociation of one of the 
colliding Au nuclei are shown in Fig. |2| as functions of the beam energy. In 
addition to Fig. |2|, where the calculations are presented only for Au + Au reaction, 
in Tab. ^| the inclusive cross sections are given for the three reactions listed in 
Tab. |||. One can see that the partial cross sections for neutron emission are 
especially large. This is due to the high average neutron multiplicities, 4.1, 7.2 
and 8.8, at SPS, RHIC and LHC energies, respectively. 

The relative contributions to the inclusive cross sections from the double pho- 
ton absorption, a?' / (o^ + cr^), are listed for the same reactions in Tab. f|. By 
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inspecting Tab. [| one can conclude that about 10% , 5% and 2% of particles, 
on average, are produced in the second order processes at SPS, RHIC and LHC 
energies, respectively. These values should be compared with the corresponding 
contributions to the total ED cross sections (Tab. ||), /<ted, which are, 3.5%, 
2.3% and 1.4%, respectively. Thus, the naive expectation, that in the double pho- 
ton absorption the multiplicities of particles should be exactly twice as large as in 
the single photon absorption, is not valid. Indeed, according to this expectation 
the partial fractions should be twice as large as the fractions of in the total 
cross section, i.e. 7% , 4.6% and 2.8% . This rule is violated for the intermediate 
mass and fission fragments. Nevertheless the general tendency is obvious: high 
order corrections become less and less important with increasing energy. This is 



true not only for the total ED cross section, as it was recently noticed in Ref. [[35 
but also for the partial cross sections for nucleon and fragment emission and pion 
production. 



3.2 Neutron multiplicity distributions 

As stressed above, neutrons are the most abundant particles produced in the ED 
process. In addition, there are many mechanisms of neutron emission that makes 
the calculation of their characteristics especially difficult. During the fast stage 
of the photonuclear reaction some intranuclear nucleons may escape from the 
nucleus via a direct knock-out either by the incoming photon or by a secondary 
particle in the course of the intranuclear cascade. In addition, as it was mentioned 
before, many neutrons may be evaporated from an excited residual nucleus. 

The neutron multiplicity distributions predicted by the RELDIS code are 
presented in Fig. |3] for the three reactions discussed above. One can see that 
these distributions have a non-trivial structure. They are strongly peaked at the 
In emission channel associated with the GR decay. On the other hand, there 
is a long tail of multiple neutron emission associated with the knock-out and 
evaporation processes. One can see, for example, that the probability to emit 
more than 20 neutrons is quite noticeable (~ 5% at RHIC). The mean neutron 
multiplicities for the three reactions considered are 4.2, 7.2 and 8.8 at SPS, RHIC 
and LHC energies, respectively. These results might be important for designing 
zero-degree calorimeters at RHIC and LHC. 



3.3 Spectra of nucleons, pions and nuclear fragments 

Now let us discuss the differential distributions in rapidity yi a b and transverse 
momentum P t for nucleons, nuclear fragments and pions. 

The Lorentz- invariant inclusive differential cross sections, da/dyi a b and da/P t dP t , 
for neutrons are shown in Fig. |4]. Obviously, the dominant part of neutrons is con- 
centrated near the beam rapidity (see upper panel). However, the cross section of 
neutron emission is quite large (~ 0.1 — 1 b) even 2-3 units away from the beam 
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rapidity. In our model these neutrons are produced due to reactions 7iV — > inn 
(1 < % < 8) initiated by a high-energy virtual photon. The recoil neutron in these 
reactions has a chance to keep a great part of the photon energy. 

The neutron P t distributions (lower panel) reveal clearly the presence of two 
components associated with fast and slow stages of the reaction. The evaporation 
neutrons are concentrated mainly at P t < 150 MeV, while the high momentum 
tail of the P t - distribution is formed by early emitted knock-out neutrons. 

Fig. |^ shows the invariant differential cross sections, da/dyi^ and da / P t dP t 
for protons. Although the proton rapidity distributions are similar to those for 
neutrons, their transverse momentum spectra reveal a difference. Namely, the 
evaporation component in the proton spectra is suppressed by the Coulomb bar- 
rier. Instead of two components, a broad peak is seen at P t ~ 150 — 250 MeV/c. 
Its origin can be understood as follows. A virtual photon from the region of 
quasideuteron absorption, E 1 ~ 30 — 100 MeV, produces a pair of nucleons with 
approximately equal momenta of 200-300 MeV/cQ. They are emitted predomi- 
nantly at angles 80° — 85° with respect to the collision axis in the nuclear rest 
frame and therefore the corresponding P t is approximately 150—250 MeV/c. This 
prediction can be checked by measuring the P t distributions of protons produced 
in the ED reactions. Up to now mainly the collective excitations of nuclei in the 
GR region were studied in these reactions through the neutron emission []27] or 
fission [[36] channels. There were no special attempts to see how virtual photons 
interact directly with individual nucleons. 

The transverse momentum distributions of intermediate-mass (3 < Z < 30) 
and fission (30 < Z < 50) fragments are shown in Fig. ^. They are remarkably 
flat. All these fragments are produced at late stages of the reaction and their 
characteristic momenta, in the nuclear rest frame, are of about 100 MeV/c per 
nucleon. Their rapidity distributions are not shown because they are extremely 
narrow and centered around the beam rapidities. As shown in Ref. ||, shell 
structure effects do not play a role at these high excitation energies and, there- 
fore, fission fragments have a symmetric charge distribution. The distribution of 
intermediate-mass fragments is close to a power law. 



3.4 Background for 77 processes 

In this section we demonstrate that pions produced in the 7 A process may cause 
problems in studying 77 fusion reactions. The predicted distributions of positive 
and negative pions from jA processes are given in Figs. [7| and |8|, respectively. 
Both 7r + and n~ have very wide rapidity distributions which extend toward mid- 
rapidity up to \yiab\ < 2 for LHC and up to yi a b ~ for RHIC. One can also notice 
that some pions have rapidities which exceed the beam rapidity {vrhic — 5.4, 
Ulhc — 8-3). This is a peculiarity of the multiple pion photoproduction at high 

1 Fermi motion is neglected in this estimation. 
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energies when a large phase space is available for pions. Some pions are emitted 
backward with respect to the virtual photon momentum. 

The pions produced in jA interactions close to the mid-rapidity region may 
create complications for the detection of particles produced from 77-fusion pro- 
cess [§, |9], [UJ. This can be seen from a simple estimate. As predicted in Ref. [|Kj, 
the cross section to produce a £)+ meson at LHC is 0.144 mb and da/dy^ is 
peaked at yi a b = so that almost all these mesons will have rapidities within the 
interval \yi a b\ < 2. The probabilities of D£ decays into r]7i + , uj-k + or p°7r + chan- 
nels are expected at the level of 1 — 5%, [[37fl. Then the cross section to produce 
a 7r + from the decay turns out to be below 0.007 mb. On the other hand, for 
7r + from the jA process we predict da/dyi^ ~ 10 mb at \yi a b\ = 2 (see Fig. |7]). 
These estimations show that it will be very difficult to separate the products of 
the .D+ decay from the background of the 7 A process. This does not necessarily 
mean that the detection of Df is impossible since different kinematical cuts may 
be applied to suppress the •yA background, as well as other decay modes may be 
used for the detection. The considered example rather demonstrates the im- 
portance of the proper evaluation of the jA background in planning experiments 
on 77 fusion at future colliders. 



4 Conclusions 

In this paper we have further developed the model of electromagnetic dissociation 
of nuclei first introduced in Ref. ||. The computer code RELDIS is developed 
to calculate a wide range of observable characteristics associated with the ED 
reactions. This code includes the dissociation channels with nucleons, pions, 
fission fragments and intermediate mass fragments from nuclear multifragmen- 
tation. The model calculations are in a good agreement with experimental data 
obtained for the dissociation of 1 GeV/nucleon 197 Au ions. This gives us the 
basis to put forward predictions for the electromagnetic dissociation of ultrarel- 
ativistic heavy ions at RHIC and LHC energies. One should bear in mind that 
the multiplicities of particle produced by high-energy virtual photons are rather 
small as compared with those in central nuclear collisions. But the reaction rates 
are large. Besides the fundamental interest, the ED process may have interesting 



applications. As proposed in Ref. Pql , the neutron emission in the ED process 
may be used for luminosity monitoring via zero degree calorimeters. Detailed 
calculations presented in this paper make it possible to eliminate some unneces- 
sary simplifications in evaluating the neutron transverse momentum distributions 



which were adopted in Ref. f38|. Although the advantages of detecting mutual 



one-neutron dissociation events are prominent [3J|, one should carefully estimate 



the contribution of more deep dissociation processes, as it has been made in the 
present paper. The electromagnetic fission of Au and Pb nuclei will be also no- 
ticeable at RHIC and LHC. An alternative approach to luminosity monitoring 
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would be to measure the yields of fission fragments or protons, which can be 
easily separated from the beam nuclei by their different Z/A-ratio. Using pre- 
sented above Pt distributions of fission and intermediate mass fragments one can 
estimate the response of forward detectors to such fragments. The calculated 
differential distributions of pions should be taken into consideration in planning 



experiments on 77 fusion at CERN [40 



One can use the predictions of the present model in order to disentangle differ- 
ent mechanisms of the heavy-ion dissociation, not only though the measurement 
of the Z-dependence of the total cross section, as it was made in Ref. |4l| , but 
also through the registration of some specific exclusive channels (multiple neu- 
tron emission, for example). Such measurements are needed, in particular, for 
understanding the origin of the correction term in the dissociation cross section 



found in Ref. [41 
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Table 1: Calculated and measured |^7| cross sections for nucleon emission in the 
reactions 197 Au( 197 Au, 



,197-rc 



Au and 197 Au( 197 Au,pxn) 196 - a; Pt at 1A GeV. 





Cross section (mb) 


Channel 


Model 


Experiment 








0i 






°P 




n 


2432 


74 


2506 


3077 ± 200 


2n 


926 


13 


939 


643 ± 105 


3n 


129 


77 


206 


171 ± 26 


Total 










n — 3n 


3487 


164 


3651 


3891 ± 227 


An 


40 


43 


83 




5n 


7.7 


11.6 


19.3 




6n 


2.2 


4.4 


6.6 




Total 










An — 6n 


49.9 


59 


108.9 




P 


1.7 





1.7 




pn 


3.2 





3.2 




p2n 


2.5 


0.6 


3.1 




Total 










p — p2n 


7.4 


0.6 


8 




p3n 


0.8 


0.7 


1.5 




pAn 


0.3 


0.5 


0.8 




p5n 


0.1 


0.2 


0.3 




Total 










p3n — p5n 


1.2 


1.4 


2.6 
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Table 2: Total cross sections (barn) of electromagnetic dissociation, (Jed, an d 
contributions from the single, and double, <j( 2 \ photon absorption. The 

values are given for the present and future heavy-ion beams of Pb and Au nuclei. 



Reaction 


a™ 




&ED 


158A GeV 

208p b on 208p b 


41.2 


1.5 


42.7 


100A+100A GeV 
197 Au on 197 Au 


93.8 


2.2 


96 


2.75A+2.75A TeV 

208p b Qn 208p b 


208.7 


3 


211.7 



Table 3: Inclusive cross sections (barn) of emission of nucleons, nuclear fragments 
and pions in the electromagnetic dissociation of Pb and Au nuclei calculated for 
the present and future RHI beams. 



Reaction 




P 


n 


3 < Z < 30 


30 < Z < 50 


7T + 


7T 




158A GeV 

208p b on 208p b 


15.9 


175 


0.09 


0.5 


1.37 


2.28 


2.7 


100A+100A GeV 
197 Au on 197 Au 


104.7 


686.7 


6.28 


4.3 


24 


30.2 


33.1 


2.75A+2.75A TeV 

208p b Qn 208p b 


302.4 


1853 


22.3 


11.9 


90.4 


111.6 


121.4 



Table 4: Electromagnetic dissociation of Pb and Au nuclei. Contributions (% ) 
due to the double photon absorption to the inclusive cross sections of emission 
of nucleons, nuclear fragments and pions. 



Reaction 




V 


n 


3 < Z < 30 


30 < Z < 50 




7T 




158A GeV 

208p b on 208p b 


8.1 


6.7 


14. 


10.5 


11.6 


11.1 


10.1 


100A+100A GeV 
197 Au on 197 Au 


3.7 


3.3 


4.2 


4. 


5.6 


5.3 


5.06 


2.75A+2.75A TeV 

208p b Qn 208p b 


1.8 


1.8 


1.4 


1.8 


2.4 


2.4 


2.3 
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Figure 1: a) Ratio of average excitation energy of residual nucleus E* to input 
photon energy E 7 , b) absolute values of E* per nucleon of residual nuclei, c) aver- 
age photoneutron multiplicities (N n ) and d) width of photoneutron multiplicity 

W n = J (N%) — (N n ) 2 as functions of E 1 in photoabsorption on lead nucleus. 
Solid lines show the RELDIS calculations, dashed lines and point represent the 



experimental data |25 
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SPS RHIC LHC 




Figure 2: RELDIS predictions for inclusive cross sections of emission of nucleons, 
nuclear fragments and pions in the electromagnetic dissociation of one of the 
colliding Au nuclei. Cross sections for intermediate mass (3 < Z < 30) and 
fission (30 < Z < 50) fragments are labelled as "IMF" and "FIS", respectively. 
Thick dotted line shows the total electromagnetic dissociation cross section, y/s 
values for SPS, RHIC and LHC are shown by arrows. 
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Figure 3: RELDIS predictions for multiplicity distributions of neutrons in the 
electromagnetic dissociation of Pb nuclei at LHC and SPS energies (solid and 
dotted lines, respectively) and Au nuclei at RHIC energies (dashed lines). 
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Figure 4: RELDIS predictions for inclusive rapidity (top panel) and transverse 
momentum (bottom panel) distributions of neutrons emitted in the electromag- 
netic dissociation of Pb nuclei at LHC and SPS energies (solid and dotted lines, 
respectively) and Au nuclei at RHIC energies (dashed lines). Rapidity distri- 
butions for colliders were reflected respect to the value yi a b = 0. Arrows show 
rapidity of corresponding beams. 
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Figure 6: RELDIS predictions for inclusive transverse momentum distributions 
of intermediate mass (top panel) and fission (bottom panel) fragments emitted 
in the electromagnetic dissociation of Pb nuclei at LHC and SPS energies (solid 
and dotted lines, respectively) and Au nuclei at RHIC energies (dashed lines). 
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